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Abstract 

We present a lattice QCD calculation with two dynamical flavors of the isovector vector 
correlator in the high-temperature phase. We analyze the correlator in terms of the 
associated spectral function, for which we review the theoretical expectations. In our 
main analysis, we perform a fit for the difference of the thermal and vacuum spectral 
functions, and we use an exact sum rule that constrains this difl'erence. We also perform 
a direct fit for the thermal spectral function, and obtain good agreement between the 
two analyses for frequencies below the two-pion threshold. Under the assumption that 
the spectral function is smooth in that region, we give an estimate of the electrical 
conductivity. 
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1 Introduction 



The properties of strongly interacting matter under extreme conditions are the subject 
of intensive experimental and theoretical investigation. A comprehensive picture of a 
state of matter requires not only the knowledge of equilibrium properties such as the 
equation of state and static susceptibilities, but also an understanding of its transport 
properties. 

In the high-temperature phase of QCD, the transport coefficients (the shear and 
bulk viscosities as well as the electrical conductivity) have been calculated perturba- 
tively to full leading order in the strong coupling fl}|3]. However, in the range of 
temperatures that can be reached in heavy ion collisions, the perturbative uncertainty 
remains large. Phenomenologically, the observation of large elliptic flow in heavy ion 
collisions at RHIC and at the LHC hints at a small shear viscosity ((4|[5] and references 
therein). Furthermore, the measured spectrum of dileptons is related to the spectral 
functions of the electromagnetic current, integrated over the history of the expanding 
system [6]. 

Any reliable calculation of a transport coefficient of QCD at a temperature of a few 
hundred MeV would be extremely valuable. It therefore makes sense to undertake a cal- 
culation in the comparatively easiest possible channel. For a lattice QCD approach, the 
spectral function of the isovector vector current is perhaps the most accessible channel. 
First, the correlation function is evaluated via a single, connected Wick contraction, 
allowing for a good signal-to-noise ratio in the Monte-Carlo simulation. Second, the 
light quarks do not introduce a new dynamical scale into the problem in the way that 
heavy quarks do, which typically requires the use of an effective field theory. Third, the 
corresponding vacuum spectral function, being extremely well known experimentally 
due to decades of measurements of the R ratio and of r decays (see for instance [?]), 
provides a useful reference. 

Here we present the first lattice calculation of the Euclidean isovector vector corre- 
lator in the high-temperature phase of QCD with dynamical quark flavors and analyze 
it in terms of the spectral function. We adopt an approach used previously in the bulk 
channel jSj, which consists in analyzing directly the diff'erence of the thermal and vac- 
uum correlators. Moreover, we exploit a recently derived sum rule which constrains 
the integral over the difference of spectral functions (divided by frequency) to vanish. 
We compare our results at finite lattice spacing with a recent analysis performed in the 
continuum limit of quenched QCD [o] . 

In spite of the technically favorable properties of the channel, determining the 
vector spectral function with frequency resolution Au <C T remains a numerically ill- 
posed problem (see for instance the discussion in [lO|). Our main goal in this paper 
is therefore of qualitative nature and consists in determining the gross features of the 
thermal spectral function, in particular in which frequency bins (of width Ao; > 2T) it 
under- or overshoots the vacuum spectral function. 
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2 Theoretical expectations for the spectral function 



In this section we set up our notation and define the relevant correlation functions. We 
summarize the theoretical expectations for these correlators and the associated transport 
properties. 



2.1 Definitions 

Our primary observables are the Euclidean vector current correlators, 

G^,{t,T) = J d^x{J^{T,x) MO)^), (2.1) 

with J^{x) = [u{x)jij,u{x) — d{x)'j^d{x)) the isospin current. The expectation values 
are taken with respect to the equilibrium density matrix e~^^ /Z((3), where (3 = 1/T is 
the inverse temperature. The quark number susceptibility is defined as 



Xs = - J <fx (Jo(x) Jo(0)) = -p j d'x ( Jo(r, a;) Jo(0)). (2.2) 

Due to charge conservation, the two correlators of interest are exactly related via 

Gu{t,T) = XsT + G^^{t,T) . (2.3) 
The Euclidean correlators have the spectral representation 



oo 



r (^T\ / n (,, T\ cosh[a;(/3/2 - r)] 

^^^^^^ I 2^ P^^^""' ^) sinh(u;/3/2) " ^^'^^ 







For a given function p{uj,T), the reconstructed correlator is defined as 

°° . ^/N cosh[w(f - r)] 



G-^(T,r;r')= / — piu:,T') . - (2-5) 

Jo 27r smh(Li;p/2) 

It can be interpreted as the Euclidean correlator that would be realized at temperature 
T if the spectral function was unchanged between temperature T and T' . For T' = it 
can be directly obtained from the zero-temperature Euclidean correlator via [S] 

G'^^{t,T) = G'^^{t,T;0)= ^ G(|T + m/3|,r = 0). (2.6) 



2.2 Theoretical predictions 

The isospin diffusion constant D is given by a Kubo formula in terms of the low- 
frequency behavior of the spectral function, 

1 pii{uj,T) 

Dxs = T lim • (2-7 

6 t^->-o a; 
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In the thermodynamic hmit, the subtracted vector spectral function obeys a sum 
rule (see [IT] sec. 3.2), 



duj 



ijj 



Ap{oj,T) = 0, 



Ap{uj,T) = pii{u,T) - pii(w,0). 



(21 



This sum rule is based on two ingredients. Firstly, the two-point function of a spatial 
component of the vector current at vanishing four-momentum can be interpreted as 
the susceptibility of the isospin charge at zero temperature in a system with one short 
spatial periodic dimension of length (3. As long as the correlation lengths are finite, this 
susceptibility vanishes. Secondly, subtracting the same quantity in the infinite-volume 
vacuum enables one to write a convergent sum rule in u for the susceptibility. It is well 
known from the operator-product expansion that the difference of spectral functions in 



12 



13 for explicit calculations. 



Eq. (2.8) falls off as l/w^ at large frequencies, see 

For non-interacting massive quarks in the fundamental representation of the SU(A'^c) 
color group, the vector spectral function is diagonal in flavor space and takes the form 



zvr 



(2.9) 
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up' tanh 



4T' 



The next-to- leading order has been computed very recently [14] . At large frequencies the 
radiative corrections (1 + Qs/vr + . . . ) to the coefficient of the uP term are temperature 

)■ 



16 



independent and known to order [15] (for quark mass effects in the vacuum, see 
The susceptibility and the mean squared transport velocity (f^) have the following 
expressions. 



4iVc/3 
4iVc/3 



S'Tp 
(2^ 



/p(l - fv) 



(27r)3 ^P^^ ^2 



(2.10) 
(2.11) 



with /p 



\l\e^^^ + 1] the Fermi-Dirac distribution and Ep = \/ p'^ + m?. 



It is now 

straightforward to check that the sum rule (2.8) is verified in the free theory. The 
positive contribution of the transport peak (the uj5{bj) contribution in the free theory) 
is compensated by a deficit of the thermal spectral function at intermediate frequencies, 
LO = 0(T) in the massless case. The susceptibility and mean square velocity have simple 
expressions in the massless and in the heavy-quark limits, 





m = 


m ^ T 


Xs 


Ncrp2 
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T 


mT 
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3T/m . 





(2.12) 
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Beyond the non-interacting theory, at weak couphng kinetic theory predicts the 
presence of a narrow transport peak in the spectral function at a; = 0, whose width 
and height are related to the properties of the quasi-particles. Introducing a separation 
scale A between the transport time scale and the thermal time-scale, the area under the 
transport peak is, to leading order, preserved by the interactions (TtI, 



AA)= ["p^-^^ = Xs{v\ (2.13) 
The width of the transport peak however becomes finite. In the heavy-quark limit for 



instance, the Langevin effective theory predicts a Lorentzian form 17 



Pii{oj,T) = Xs{v'')^P^, m:^T, (2.14) 

where r] is the 'drag coefficient', 1/^/ = D^- The case of massless quarks can be treated 
with the Boltzmann equation [T| |18|19] , with a form of the spectral function qualitatively 



similar to Eq. (2.14) and a drag coefficient given by 



r]=^CFmllog{T/mD), (2.15) 

with mD the Debye mass and Cp = 2N ■ 

Finally, in contrast with the weak-coupling analysis outlined above, it is worth 
mentioning that at least one theory is known where the vector spectral function does 
not exhibit a transport peak. In the strongly coupled = 4 super- Yang-Mills theory, 
the spectral function of the R-charge correlator reads |^| , 

Ri rj.. 3xs sinh^ 
fe(-,T) = — [-) ^^^tr^_^^^^ . (2.16) 

The static susceptibility is given by Xs = —% — , and the diffusion constant hy D = 



2.3 Connection with electromagnetic observables 

The electrical conductivity a is extracted from the correlator of the current J^™ = 
'Ylif QfQflfiQf- If the quark species are degenerate, this correlator is given by the sum of 
a contribution proportional to and a contribution proportional to (^j Q/)^- At 

high frequencies, the latter contribution to the spectral function is predicted to be small 
in perturbative QCD. Assuming that this contribution is also small at low frequencies, 
we have for the electrical conductivity in A^f = 2 QCD 

a = CeraDXs- (2.17) 

with Cem = 'Ylf=udQ'^f = 5/9. If one further assumes that a is fairly insensitive to 
the quark mass, and that the transport properties are not significantly affected by the 
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presence of virtual strange quark pairs, the electrical conductivity in A'^^ = 2 + 1 QCD 



is given by Eq. (|2.17|) with Cem = Ef=uds = 2/3- 



f=u,d,s ^/ 

ogical significance of the electromagnetic spectral function is that 



The phenomeno 

the dilepton rate and the real-photon production rate are given by |18l|21 



dud^p 67r3 (a;2 _ p2) (g^/T _ i) ' 

dR 3 

lima;-^ = —.a{T)Taem (2.19) 
w^o d^p 27r2 

where Oem is the electromagnetic fine structure constant. 
2.4 Phenomenology of the p resonance 

In the vacuum, the QCD spectral function of the electromagnetic current is well mea- 
sured via the R{s = cj^) ratio. The p meson completely dominates the spectral function 
up to about uj = IGeV. We work in the exact isospin symmetric theory and therefore 



ignore issues related to isospin breaking and p — u) mixing, see 22 for a recent reference. 

Since we are working with the isospin current, we should restrict ourselves to final 
hadronic states with 1 = 1. The Ri{s) ratio is defined analogously to R{s) with this 
restriction. A rough parametrization of the experimentally measured R\{s) ratio was 
given in Ref . |11| , Eq. (93) . The spectral function in our normalization is related to the 
i?i ratio via 

Pii{oj,Q) = -Ri{u:^)uj^. (2.20) 
vr 

We can now make a simple argument about the thermal spectral function based on the 
exact sum rule (2.8), the kinetic theory sum rule ( 2.13[ ) and the experimentally known 



vacuum spectral function. Using the parametrization of |lT], the area under the vacuum 
spectral function up to ci; = IGeV is about 

r'"-^^»CH4GeV^ (2.21) 



.,0 vr w 

By contrast, the corresponding area for free massless quarks at zero temperature is 

Aree = O.OTGGeV^. (2.22) 

Taking into account the physical light quark masses changes this value by a negligeable 
amount. 



In the free theory, the sum rule (2.8) is satisfied. Now switching on interactions 
between quarks, let us assume for the sake of the argument that in the high temperature 
phase they can be described perturbatively. In the small frequency region, interactions 



turn the delta function in Eq. (2.9) into an approximate Lorentzian curve |18j, preserving 



its area to leading order. This is the content of a kinetic theory sum rule 17 . In the 



vacuum, interactions have a dramatic effect on the spectral function, due to chiral 
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Figure 1: The phenomenological isovector vector spectral function pa in the vacuum compared to 
the leading-order weak- coupling prediction for the spectral function in the high-temperature phase. 



The Born term is as given in Eq. (2.9) for m = and we have represented the transport peak 
around uj — by a Lorentzian with a width matched to fl^ for as — 0.3 and an area determined 
by {v'^) = 1 together with Xs/xt^ — 0-88 and 0.9 respectively for the left and right panel (table^^ 
and \23^ ). The thin red line represents the vacuum spectral function for non-interacting massless 
quarks. The weak-coupling thermal spectral functions must receive additional contributions of 



the size given by the 'missing area' rectangle in order to satisfy the sum rule (2.8) 



symmetry breaking and confinement, and convert its area up to IGeV from ^free to An. 



Since the sum rule (2.8) must still be satisfied, the weak-coupling spectral function in 



the high-temperature phase must acquire an additional area of 

-4missing = Ap- Arcc ~ O.OSSGeV^. (2.23) 



Note that at very high temperatures, the area (2.23) is negligible compared to the area 
^(A) = Xs{v'^) under the transport peak, which grows as T^. However, consider the 
situation at T = 253MeV (the temperature at which we have computed the correlator 
on the lattice, see the next section). Assuming the existence of a transport peak, its area 
^(A) is about O.OSeGeV^, if we correct for the fact that the static isospin susceptibility 
is about 12% below its Stefan-Boltzmann limit x^b (^ggg table [l] and [23]). The area 
missing from the weak-coupling spectral function is thus comparable in size to the 
transport peak area at this temperature. The argument is illustrated in Fig. ([T]). 

At cj ^ T, the difference of spectral functions Ap{uj, T) can be analyzed with 



the operator product expansion (see (13, 24 and References therein). The lowest' 



dimensional gauge- invariant operators are of dimension four, and therefore the asymp- 
totic behavior is Ap{uj,T) oc l/w^ (possibly up to logarithms). According to Eq. (4.1) 



of Ref. 24 , the leading term of order is positive. However, at T = 253MeV its 



contribution to the area under Ap{uj,T) /oj is too small to explain the missing area 



( 2.23D . 

In conclusion, at temperatures that are accessible in heavy-ion collisions, the sum 
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rule (2.8) and the R{s) ratio measurements place an important constraint on the thermal 
spectral function. 



3 Lattice QCD data 

All our numerical results were computed on dynamical gauge configurations with two 
light, mass-degenerate 0(a)-improved Wilson quark flavors. The configurations were 



generated using the MP-HMC algorithm 25 , 261 in the implementation of Marinkovic 



and Schaefer [27] based on Liischer's DD-HMC package [28] . The zero-temperature 
ensemble was made available to us through the CLS effort 29 . The gauge action is 
the standard Wilson plaquette action |30| , while the fermions were implemented via 
the Wilson-Clover discretization with non-perturbatively determined clover coefficient 



31 



We calculated correlation functions using the same discretization and masses as 
in the sea sector in two different ensembles. The first corresponds to virtually zero- 
temperature on a 64^ x 128 lattice (labeled 07 in [32|) with a lattice spacing of a = 



0.0486(4) (5)fm [32] and a pion mass of = 270MeV, so that m^L = 4.2. Secondly 
we generated an ensemble on a lattice of size 64^ x 16 with all bare parameters iden- 
tical to the zero-temperature ensemble. In this way it is straightforward to compare 
the correlation functions respectively in the confined and deconfined phases of QCD. 
Choosing Nr = 16 yields a temperature of T ~ 250MeV. Based on preliminary re- 
sults on the pseudo-critical temperature Tc of the crossover from the hadronic to the 
high-temperature phase [33], the temperature can also be expressed as T/Tc « 1.2. 
The renormalization of the vector correlator assumes the form 

G'^.(t,/3,5o) = ^y(9o)[l + bv{go)am,f (g'^^{t, ^, go) + ...). (3.1) 

The dots refer to 0(a) contributions from the improvement term proportional to the 



derivative of the antisymmetric tensor operator [34 ,35 that we did not compute. We 
note however that perturbatively, the 0(o) contribution is suppressed by two powers 
of Us- We use the non-perturbative value of Zy provided by [36]. From the perturbative 



results of 35 we estimated the magnitude of the term ~ 6y(g'o)a?Wg to be of the order 



"q 

of 0.3%. Considering that the error of Zy is at the 1% level, we drop this contribution 
altogether in the following. Note that since we use the same bare parameters at zero and 
finite temperature, the multiplicative renormalization of the vector current only affects 
the corresponding correlation functions - and their difference - by an overall factor. 

The vacuum correlator serves as a reference in this work. To fix the parameters 
of the lightest vector state in the finite volume of the simulation, we fitted the vacuum 
correlation function to an Ansatz of the form 

G,,(r, 0) = Aie-""'^ + exp(-J]r) {n^r + 2n/T^ + , (3.2) 
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Figure 2: The vacuum vector correlator Guij^T = 0) computed on a lattice sized 64'^ x 
128 and = 270MeV, labeled 07 in f32'J. The blue line denotes the fit result using 



Eq. 3.2. The insertion shows the corresponding effective mass for < r/a < 32. 



which is the Laplace transform of 

= ^1 - mi) + ^e{u: - n) Kou'. (3.3) 
The finite-volume spectral function consists of Dirac 5 functions. However, since many 



states contribute at short distances, the continuum in Eq. (3.2) should be interpreted 
as resulting from the contributions of many energy eigenstates. In Fig. [2] we show 
the vacuum vector correlator Gjj(r, T = 0) and the fit result. The insert shows the 
corresponding effective mass. The fit (performed in the interval 5 < r/a < 64), clearly 
provides a good description of the data. For r/a > 32 the quality of the data deteriorates 
and the signal is lost. 

We will only use the result for the parameters Ai and mi in the following. The 
question arises of the relation of these parameters to the infinite-volume spectral func- 



tion. An answer is given by Refs. 37,38 . Assuming phenomenologically reasonable 



values of the p coupling to the vrvr channel, the mass mi in fact appears to be within 
10% of the (infinite- volume) p mass [37] . 

In addition we estimate the thermal (isovector) quark number susceptibility Xs/T'^ 



from the time-time component of the vector correlation function (see Eq. 2.2 ), by reading 
off the midpoint of the time-time correlator. We checked the obtained value by fitting 
the correlator to a constant in the region 2 < r/a < 14 and found only a negligible 
deviation. The parameters used in our lattice setup, the 'p-meson' parameters and the 
value of the static susceptibility are summarized in Tab. [T] 
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64^ X 128 


64^ X 16 


Ref. 




5.50 






K 


0.13671 








1.7515 






Zv 


0.768(5) 


36 




a[fm] 


0.0486(4)(5) 


32 




m^[MeV] 


270 


32 




T [MeV] 
mi/T 

Kq 

n/T 


4.42(31) 
3.33(5) 
1.244(5) 
5.98(11) 


253(4) 
0.871(1) 





Table 1: The top part of the table shows the common quantities characterizing the zero- 
temperature and finite-temperature ensembles. In the lower part, the fit parameters for 
the vacuum correlator in units of T = 253MeV and the value of the (isospin) quark 
number susceptibility Xs/T"^ o-re given. For more details on the generation of the Nr = 
128 ensemble, see fsW. The number of configurations generated with N-^ = 16 is 317. 



3.1 Thermal and vacuum correlators 

In Fig. |3] we show the correlator Gii{T, T) computed at T ~ 250MeV together with the 
corresponding free 'continuum' and free 'lattice discretized' correlation functions. In 



addition we show the reconstructed correlator G[f'^(r) as obtained from Eq. 2.6 The 
results for the thermal and reconstructed correlators for r/a > 4 can be found inTab.[2j 
The reconstructed correlator lies somewhat lower than the thermal correlator. The 
insert in Fig. [s] displays the ratio Gjj(T)/G[/^(r) in order to make their relative r 
dependence visible. For small Euclidean times r < /3/4 this ratio is unity, above it 
increases monotonically until it levels off around the midpoint at about 10% above 
unity. A thermal modification of the spectral function has thus taken place (recall 
that the spectral function underlying G^f^ir) contains the bound states of the confined 
theory) . 

In Fig. |4]we show the difference 

AG(r,T) . 0..(r,T) - aTir,T) = f Ap(.. r)S2!M;f^ (3.4) 

Jo 27r smh(a;p/2j 

of the thermal and the reconstructed correlators. Given that Pn{u}, T) and Pii{uj, T = 0) 
have the same ~ behavior, this means we are subtracting non-perturbatively the 
ultraviolet tail of the spectral function. Using this difference we are therefore able to 
probe the change in the vector spectral function from the confined to the deconfined 
phase for frequencies co < 0(T). 
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Figure 3: Thermal Gii{T)/T^ and reconstructed GI['^{t)/T'^ vector correlators over Eu- 
clidean time separation r compared to the free (continuum and lattice) cases. The re- 
constructed correlator was computed by applying Eq. \2.5\ to the data obtained from a 
lattice sized Nu = 64 and Nr = 128 . The insertion shows the ratio Gii{T)/G1f^{T). 



For small times, the difference (3.4) turns out to be negative, while it is positive 
for r > /3/4. Note the errors decrease with increasing Euclidean time throughout the 
available range. We show a more detailed view of the region r/a > 3 in the insert of 
Fig. |4] Here the difference still exhibits a mild increase and levels off near the midpoint. 
The value it reaches at the midpoint is AG(r = I3/2,T)/T^ = 0.291(55). 



3.2 Comparison with the free-quark correlator 

In the previous subsection we compared the thermal correlator to its zero-temperature 
analogue. Now we analyze the thermal correlator in relation to the non-interacting 
case, which by asymptotic freedom corresponds to the regime of asymptotically high 
temperature. However, before discussing the departure of the simulation data from the 
non-interacting case, we address briefly the issue of cutoff effects. 

Since we only have data at one lattice spacing, the value of the lattice correlator, 
viewed as an estimator of the continuum correlator, is necessarily ambiguous at some 
level due to cutoff effects. To estimate the size of this ambiguity, we show the ratio 
Gu/gII^^''^^^ in Fig. 5 left), where c and I denote the analytically known free continuum 
and free lattice cases, respectively 39 . Concentrating on the ratio to the free continuum 
correlator we observe a decreasing trend throughout the entire available Euclidean time 
range, a behavior very similar to that seen in quenched studies |9j. Taking the ratio to 
the free lattice correlation function on the other hand the results are much flatter and 
almost constant at small times r/a < 3, while for r/a > 5 the two ratios track each 
other and are separated by a small shift of roughly 2%. The difference between the two 
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Figure 4: Difference AG{t,T)/T'^ of the thermal vector correlator and the correspond- 
ing reconstructed correlator as a function of Euclidean time r. The insertion shows 
AG{t,T)/T^ in the region r/a > 3. 



4 
5 
6 
7 



12.534(94) 

6.775(65) 

4.510(51) 

3.569(45) 

3.300(44) 



12.456(59) 

6.602(41) 

4.264(35) 

3.282(33) 

3.008(34) 



Table 2: Results of the thermal and reconstructed correlation functions for r/a > 4. 
Note that all results have been renormalized using the value of Zy in Table ([I]) and 
normalized by T^ . 



curves comes from the fact that the free lattice correlation function takes into account 
the tree- level lattice artifacts. 

To put this systematic uncertainty into perspective, we examine the ratio Ga/Gfi^i^^'^^'' 
in the right panel of Fig. [sj The only difference between (^^^^^'^^Z' and Gf^^f^^'^^'' is that 
the (^-function in the spectral function at w = is absent in the latter. The difference 
between the left and right panel curves thus corresponds to the contribution of the 
transport peak. It is clear from the figure that this contribution is much larger than 
the cutoff effects present at tree-level for r/a > 5. 

Returning now to the question of how much the thermal correlator differs from its 
non- interacting counterpart, we see that the simulation data lies 8 to 10% above the 
free lattice correlator. The sign of the effect corroborates our finding in section 2.4 that 
spectral weight is missing from the weak coupling spectral function. 
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Figure 5: Left: The vector correlator Gii{T) atT 250MeV normalized by the free con- 
tinuum and free discretized correlation functions Q^^^^''^^^ . Right: The vector correlator 
Gii{T) normalized by the free (continuum /discretized) correlation functions Cj/f^^'^^^ . 

3.3 Thermal moments of the correlator 

Finally we compute also the ratio of thermal moments of the correlator [9]: 

^im = 9!^ = ,Id-P(-m-^^) ^here Ki^,T) = ' . (3.5) 

Jduu'^p{uj)K{u,T) ^ ' ' 27rsinh(L^/2T) ^ ^ 

This quantity can be extracted from the correlator data by combining the results of 

iGu/G!r"'^% = /3/2) = G^/G^^^'- and: 

A,,(t)^ ^^"^^^ ^^t. = [1 + 0((/3/2 - r)^)1 ■ (3.6) 



In the free case a straightforward computation 40 yields R^^'^l = 18.423. Using this 
result together with our lattice data we obtain: 

^(0)/^(0),/.ee ^ G^^/G^'^- = 1.198(8) ^ ^ = 1.089(30) , (3.7) 

^free 



where we neglected higher-order terms in the square bracket of Eq. mWl We see that 
the ratio R^"^'^^ of thermal moments is roughly 6 — 12% larger than the free result. As 
the ratio of second thermal moments is sensitive to changes in the low frequency region 
of the spectral function (see e.g. [9]), this observation could be due to a broadening of 
the (5-function form of the free theory. 



^Including the leading correction into a fit we find it to be poorly determined by the data, while the 
constant contribution remained unchanged within errors. 
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Figure 6: The effective area ^(A, r)off, defined in Eq. {4-3), normalized by the static suscepti- 
bility. The corresponding quantity in the strongly coupled SYM theory calculated from Eq. (2.16) 
is also displayed for comparison. 



4 Analysis of lattice correlators in terms of spectral functions 

We begin with a simple but instructive analysis of the spectral function difference Ap. 



Section 4.2 contains the main analysis based on fits to the thermal part of the Eu- 



clidean correlator, and section 4.3 describes a fit directly to the thermal correlator. The 
results of the two fits are compared against each other and against previous quenched 
calculations in section 14.41 



4.1 A simple spectral analysis of the thermal part of the vector correlator 



Based on the analysis of section 2.4, it is interesting to ask whether the Euclidean 



correlator can be described and the sum rule (2.8) satisfied solely by the transport peak 
and the p-meson contribution. For this purpose, we consider the following caricature of 
Ap(a;,T), where the sum rule has already been enforced, 



Ap(a;,r) 



2tt 



C u) {6{u}) — \5{oj — mi)) , 



corresponding to the Euclidean correlator 



AG(r,T) = C 



1 mi coshmi(/3/2 — 
~j3 ~ ~2 sinhmi/3/2 



(4.1) 



(4.2) 



The mass mi is set equal to the value obtained by fitting the Ansatz (3.2 ) to the vacuum 
correlator and given in table [l] The kinetic theory sum rule (2.13) implies that C is an 
estimator for ^(A), and from Eq. (4.2) an effective value ^(A, r)efr can thus be defined 
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for every value of r, 

..A ^ AG(r,r) ,^ , 

^(A,r)cff = -, , , :v (4.3) 

\fi 2 sinhmi/3/2 j 

The result is shown in Fig. Q. The quantity ^(A, r)eff is well compatible with a 



constant value, implying that the Ansatz (4.1) already provides a good description of 



the lattice correlator. Moreover, the weak-coupling expectation that ^(A, r)efj should 
be given by Xs (up to quark mass effects reducing their average thermal velocity (f^)) 
is also well reproduced. 

It should be noted that in this simple picture the prefactor oi biyi — mi) in Eq. 



( [4l| is about -Xsmi/2 = -1A5{2)T^, to be compared with the area Ai = +4.42(31)r3 
obtained from the vacuum correlator (table [T]). This observation means that the area 
under the thermal spectral function in the region T < oj < 4T cannot be negligible 



compared to Ai, confirming the conclusions drawn from phenomenology in section (2.4). 



It is interesting to confront the lattice data with the spectral function in the strongly 



coupled SYM theory, Eq. (2.16). We therefore also plot the corresponding SYM function 
in Fig. It lies lower than the lattice data, in spite of the fact that Pj^{u}, T) / {xs^"^) 
has the same large-frequency limit as /9^p'^(w, T)/(x^''°°a;^) in QCD and also satisfies 



the sum rule (2.8). This comparison shows that the lattice data is not simultaneously 
compatible with the combination of (a) the functional form of the SYM spectral function 
and (b) the corresponding very low diffusion constant. The data is however perfectly 



compatible with the substitution of the delta function in Eq. (4.1) by a flat (or even 
oc (1 + ^w^/r^) as in the SYM case) behavior of Ap{uj,T)/u! up to about u) 4T, 
provided its area is adjusted appropriately. 

In order to parametrize Ap{uj, T) systematically, including the expected contribu- 
tions at high frequencies, we resort to the more sophisticated fits described in the next 
subsection. To anticipate the results, similar qualitative conclusions on the distribution 
of the spectral weight in A/9(a;,T) will be obtained. 

4.2 Fit to the thermal part of the vector correlator 

We proceed to investigate the behavior of the thermal part of the spectral function Ap 



by fitting the difference of the thermal and the reconstructed correlator, see Eq. (3.4). 
We know from the operator-product expansion that the difference of spectral functions 
falls off rapidly (as uj""^) for oj ^ T, and therefore focus on the region u < 0(r) in 
order to choose a fit Ansatz. As described in the previous subsection, the fact that the 
data (displayed in Fig. [4]) is positive at long distances and negative at short distances 
suggests that the thermal spectral weight exceeds the vacuum spectral weight at low 
frequencies and falls short of it at higher frequencie^ 

■^A continuum extrapolation is really needed to confirm the beliavior at short distances. 
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tanh(a)) ■ 
AG{T,T) I 







Ap(a))/a)T 


CO 

tanh(to) 









Figure 7: Left: Fits to AG{t,T)/T^ = [Gu{t) - G';^''iT)]/T^ . The blue and red results differ 
by the form of transport peak in the Ansatz. The error bands are computed from the covariance 
matrix of the fit. Right: The resulting spectral functions for both Ansdtze. 



We thus parametrize Ap using the following Ansatz for a; > 0: 
Ap{uj,T) 

PT,i(w,r) 

Apf{uj,T) 



Pt{^,T) - pb{^^,T) + ApFiuj,T), 
2cb 93 tanh{uj/T)^ 

4c w 



ioj/g)^ + V 
Pf{uj,T) - pf(w,0), 



PT,2(w,r) 

Pf{^,T) = 



4crtanh(a;/r) 

3 2 1 / ^ 

— KOJ tanh — - 
2tt \4T 



(4.4) 
(4.5) 

(4.6) 
(4.7) 



The bound state (B) and the transport peak (T) are represented by Breit-Wigner forms. 
Even such a simple Ansatz requires three parameters {cb, gs^'mB) to determine the 
bound state peak, two parameters (c, g) for the transport peak and one (k) for the 
'perturbative' contribution (F). We will therefore fix some of them using the vacuum 
correlator. In the following we set rriB equal to mi, given in Tab.[T| which we obtained 
from the exponential fit to the vacuum correlator. Note that the area under the bound 
state peak J dw pb/oo does not depend on the width in the limit where it is small. 
The sensitivity of the Euclidean correlator to the latter parameter is very small. We 
therefore perform fits for three fixed values of this parameter, and check the sensitivity of 
the result. We choose the values qb/T = 0.1, 0.5 and 1.0, corresponding to qb — 25, 125 
and 250MeV. 

The tail ~ T/a; of the Ansatz px^i violates the OPE prediction that Ap ~ (T/w)^ 
at large frequencies. It has been argued in [41| that this might lead to an overestimate 
of the transport contribution. To avoid this problem we introduce the Ansatz 2, where 
u) — )• Ttanh(a;/T). This Ansatz possesses the correct asymptotic behavior, as well as 
the expected linear behavior in lo at small frequencies. 
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Figure 8: Left panel: the impact of changing the width of the peak representing the p resonance, 
gs — 0.5T and gs — LOT. Right panel: the separate contributions pt{uj), psi^^), Pf{^) o,nd 
the full result of the fitted spectral functions for the transport Ansatz Pt,2- 



Finally, to complete the parametrization of Ap{u}), we include a weak-coupling 



term inspired by Eq. (2.9) describing the subtraction of the large frequency parts of 
the thermal and vacuum spectral functions. This contribution pp(uj,k) — )• vanishes 
exponentially as the frequency increases. 

In the next step we fit the combined Ansatze of Ap{uj,CB, gB,iT^B,c, g, n) to the 



data, while at the same time satisfying the sum rule of Eq. 2.8 to an accuracy of 
10~^. We limit ourselves to fitting the region 5 < r/a < 8 only, in order to minimize 



the influence of cut-off effects, as discussed in Sec. |3.1[ With uib determined by the 
vacuum correlator, we set qb successively to the three different values mentioned above 
and fixed k around unity, and fitted the parameters c, g and cb- The errors and error 
bands shown in the following have been computed using the covariance matrix of the 
corresponding fit for fixed values of gB and k. 

The resulting correlators and spectral functions are displayed in Fig. [7] for gB/T = 
0.50 and k = 1.10 while the fitted parameters are given in Tab. [3j In the left panel of 
Fig. [Tjthe data AG(r, T) is compared to the fits using pT,i{io) and pt,2{^) as transport 
contribution. We achieve a quasi-perfect description of the data for r/a > 4. The right 
panel shows that both Ansatze exhibit a substantial spectral weight around the origin 
and a negative contribution from the region of the p mass. 

Focussing on Ansatz 2, the left panel of Fig. [8] shows the sensitivity of the fit result 
to varying the parameter gs- Varying the width only has a small effect on the overall 
result. We also found little sensitivity to ±15% variations in k. In order to understand 
what drives the parameters to their final fitted values, we also plot separately the three 



contributions appearing in Eq. (4.4-4.7), including their respective error bands, in the 
right panel of Fig. [sj The contribution pf{^) mainly affects the intermediate frequency 
region around uj/T ~ 2. Its tail between 4<cj/T<10is largely compensated by the 
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Figure 9: Left panel: fit results for Gii{T)/Gj/'^'^{T) using both Ansdtze for the transport peak. 
The resulting ratio of thermal moments _R(-^^o)yjj(2,o)/re<, displayed on the right side of the 
plot. Right panel: The corresponding spectral functions normalized by ujT . 



tail of the Lorentzian centered at the origin, and this might well be what drives the 
width of the Lorentzian. 



4.3 Weak-coupling inspired fit to the thermal vector 

In contrast to the previous section, here we study directly the thermal vector correlator 
and its ratio of thermal moments R^'^'^\ We perform a fit inspired by the weak coupling 
form of the thermal spectral function, 

p(w,r) = pt{u:,T) + pf{u,T) , (4.8) 



where the form of the two contributions is defined in Eq. (4.6) and (4.7). At a given 
temperature this Ansatz is characterized by three parameters {c,g,K). We fit the full 
Ansatz p{uj, c, g, k) to the thermal correlator Gh^t), while at the same time demanding 
that i?^^''') be reproduced. In this analysis the three parameters c,g and k are fitted, 
and the fit range is 5 < r/a < 8 as before. 

The resulting correlators and spectral functions are shown in Fig. [9] with their fit 
parameters listed in Tab. [sj The ratio Gjj(T)/G|[^'^(r) in the left panel of Fig. [ojis well 
described by both versions pT,i and pT,2 of the transport contribution for r/a > 5, while 
also the ratio of thermal moments (given on the far right of the plot) is reproduced. For 
r/a < 5 our Ansatz fails to reproduce these points, which we suspect is partly due to 
cutoff effects. In the future it would be interesting to repeat the calculation at several 
smaller lattice spacings while keeping the temperature fixed, as in [9]. 

On the right hand side of Fig. |9]we show the resulting spectral functions divided by 
ojT. Clearly both Ansatze give very similar results that lie within errors of each other. 
The thermal correlator is even less sensitive to the asymptotic behavior of the transport 
contribution in the Ansatz than in the difference of correlators studied in Sec. 14. 2[ 
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c,g,K) 


c/T 


9/T 




>l(A = 1.5i j /I 


/„,2\ 


PT,i(w 


c, g) 
c, 5) 


0.61(10) 
0.59(16) 


1.22(26) 
5.7(10) 


1.42(188) 
2.92(228) 


0.702(201) 
0.764(244) 


0.806(231) 
0.877(280) 


piuj,c, 




c/T 






^(A = i.5r)/r2 




PT,l{uJ 


c, 5) 
c, 5) 


0.75(4) 
0.74(5) 


0.71(9) 
0.98(13) 


1.186(27) 
1.192(26) 


0.818(50) 
0.858(56) 


0.939(57) 
0.985(64) 



Table 3: Parameters obtained from both studies. Top: fitting Ap(uj, c, g, k) to AG(r, T) with 
gs = 0.5T and k — 1.10. Bottom: fitting Gh{t) to p{uj, c, g, n). In both eases fits were done 
for two different Breit-Wigner peaks in the low frequency region, pT,i{uj) ~ w and pT,2{^) ^ 
Ttanh(a;/T). Additionally the resulting area of the transport region A(A) and the mean squared 
velocity {v^) is given. For details see the text. 



4.4 Discussion 



We now compare the results of the fits to AG(r, T) and Gii{T,T). Since the vacuum 
spectral function vanishes below 2m,r ~ 540MeV (in infinite volume), pii{uj,T) and 
Ap{u},T) should be equal for w < 2m7r ~ 2. IT. We thus plot the spectral functions 



obtained from the two fits in this frequency region, see Fig. 10 Here we restrict ourselves 
to showing these results based on pT^2{^,T), as their theoretical foundation is more 
sound than those with pt,i{uj,T). All curves are multiplied by a factor 1/6, which 
means that the intercept at a; = yields an estimate of a/CemT with a the electrical 
conductivity of the quark gluon plasma. 

The results obtained by fitting Gii{T,T) agree very well with the central values 
obtained by fitting AG{t,T), as summarized in Tab. [3j The fit to AG(r, T) using the 
transport Ansatz pT,2i^) yields a slightly lower intercept. If we assume the spectral 



function to be as smooth around the origin as Fig. 10 suggests, we obtain the following 
estimate for the electrical conductivity of the quark gluon plasma at T ~ 250MeV, 

77^ = 0.40(12), (4.9) 

where Cem = Ef=u4 Q) = 5/9 for A^/ = 2 and C^m = Ef=uAs Q) = 2/3 for Nf = 2+1, 
see Sec. |2.3[ It should be remembered that the Euclidean correlator can be perfectly well 
described by an infinitely narrow transport peak (corresponding to an infinite electrical 
conductivity), see section 4.1 Although obtained under a strong assumption, it is 
interesting to compare (4.9) to other lattice results obtained under similar assumptions. 
The following comparison is made with quenched results, since to our knowledge there 
are no previous dynamical QCD studies. 

A quenched calculation using staggered fermions based on fitting the Fourier trans- 
form of the correlator obtained a/T = 7Cem ^42^1 at 1-5 < T/Tc < 3.0, where the pure 
SU(3) gauge theory critical temperature is around 290 MeV. A further study using 
staggered fermions and an analysis based on the maximum entropy method obtained 
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Figure 10: Comparison of the spectral functions obtained from analyzing (a) AG(t, T) and (h) 
Gii{T,T) using in both cases pT^^i^^) ^ tanh(cj/r) in the low frequency region. All curves have 
been multiplied by a factor 1/6 and divided by ujT , entailing that the intercept at uj = Q yields 
an estimate of a/C'emT. 



a/T = 0.4(1) Cem [43]. Finally, a recent quenched study using Wilson- Clover fermions 
in the continuum limit obtained 0.33Cem < (t/T < ICem at T ~ 1.45Tc [o]. Our results 
using dynamical Wilson-Clover fermions at Nr = 16 are thus completely compatible 
with the recent quenched results. 

Even though the transport contribution extracted from our fits is not narrow and 
an interpretation in terms of kinetic theory (Eq. 2.13) is no longer rigorously motivated, 
we also compute the effective mean squared velocity of the quarks (v^). Choosing the 
scale parameter A/T = 1.5 and numerically integrating the fitted spectral functions we 
obtain ^(A = 1.5T). Assuming Eq. (2.13) and using the quark number susceptibility 
as given in Tab. [T| it is straightforward to estimate an effective mean squared velocity 



{v )cs- The results for ^(A = 1.5T) and (t" )eff are listed in Tab. [3j For all fits we 
obtain reasonable values for )efr in the range 0.80 < {v'^)es < 0.99. Thus, although 
we are not able to demonstrate or exclude the validity of the kinetic theory description, 
the values of the effective quark velocity extracted in this way are in line with its 
expectations. By contrast, the AdS/CFT spectral function (2.16) (which clearly cannot 
be described by kinetic theory) yields an effective quark velocity of {v'^)cS = 0.47 for 
A = 1.5T. 

Using the result for pii{uj,T) from the fit to Gii{T,T), it is straightforward to 
compute the production rate of thermal lepton pairs in the quark gluon plasma with 
two light dynamical quark flavors from Eq. (2.18). The resulting rates are shown in 
Fig. 
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We give the low frequency behavior of dNi+i- /dud'^p obtained from pii{uj,T) 
using the Ansatz 2 in the transport region. Comparing our results with the free (Born) 
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1e-06 




Figure 11: The production rate of lepton pairs from two-flavor lattice QCD for frequencies lo 



given in units [GeV] as calculated from spectral functions of the vector current using Eq. 2.18 
The curve comes from fitting Pii{uj, T) to Gii{T, T) using the transport peak Ansatz 2. The black 
line shows the result from HTL perturbation theory with mr/T = 1, while the dashed line denotes 
the (free) Born rate. 



rate (dashed line) and the hard thermal loop result (black line) with a thermal mass of 
uit/T = 1 j44j, we observe that for frequencies u/T < 1.5 our result is below the result 
from HTL, above this value however it follows it very closely. 

The results shown in Fig. [TT] correspond to a system at thermal equilibrium with 
T ~ 250MeV. To make contact with results obtained in heavy ion collisions one has to 
take into account the real-time evolution of the volume of the system, using a hydro- 
dynamic model along the lines of [45| (for a recent study of out-of-equilibrium photon 
and dilepton production using the AdS/CFT correspondence see (46|[47| ). 

5 Conclusion 

In this paper we have obtained the isovector vector correlator in the high-temperature 
phase of two-flavor QCD at T ~ 250MeV as well as in the vacuum at the same set 
of bare parameters. This allowed us to analyze both the difference of the thermal 
and the vacuum correlator and the thermal correlator directly. In the former case the 
analysis is further constrained by an exact sum rule. Given the uncertainties inherent 
in trying to extract information on the spectral function from Euclidean correlators, 
the two methods give a consistent picture of the thermal spectral function in the low to 
moderate frequency range uj < 1.5T. 

The vacuum spectral function is known to receive a very large contribution from the 
p meson from experimental e~^e~ and r decay data. The main qualitative lesson we have 
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learnt is that the reduction or complete absence of such a peak and the appearance of a 
substantial spectral weight in the low-frequency region provide a very good description of 
the lattice data and are compatible with the sum rule. Moreover the area under the latter 
spectral weight matches the expectation of kinetic theory. This picture is corroborated 
by a simple phenomenological study, presented in section 2.4 based on the experimental 
R{s) ratio and the sum rule. We also note that the analytic result (2.16) in the strongly 
coupled limit of = 4 super- Yang-Mills theory exhibits a qualitatively similar (but 
quantitatively different in amplitude) change of sign in the difference of thermal and 
vacuum spectral functions, even though the theory is conformal and therefore exhibits 
no analogue of the p meson. A similar pattern was also observed in the bulk channel of 
the pure SU(3) gauge theory 

The analysis presented in this paper is based on data at finite lattice spacing. 
Obviously the next step would be to repeat the analysis in the continuum limit, as has 
been done in the quenched theory [48] . In this respect the results obtained here should 
be regarded as preliminary. We note however that our results are quantitatively quite 
close to those obtained in jis]. 

It would be very interesting to repeat the analysis carried out here in a temperature 
scan through the smooth phase transition. This would allow one to track the fate of the 
p meson from the low-temperature to the high-temperature phase and perhaps to shed 
light on the excess of dileptons observed by the PHENIX collaboration in relativistic 
gold-gold collisions [49) . The methods employed here to constrain thermal spectral 



functions may be useful in the context of cosmology as well, see for instance 50 
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